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Abstract 

We present numerical evidences for the existence of rotating black hole solutions in d-dimensional 
Einstein-Maxwell theory with a cosmological constant and for d odd. The metric used possesses 
(d+ l)/2 Killing vectors and the solutions have (d — l)/2 equal angular momenta. A Schwarschild- 
type coordinate is used for the radial variable and both signs of the cosmological constant are 
emphasized. Several properties of the solutions are studied, namely their surface gravity, mass and 
angular momentum as functions of two parameters : the magnetic field and the angular velocity 
at the horizon. The influence of the electromagnetic field on the domain of existence of the black 
holes is studied are compared to the vacuum case where analytic solutions are available. 
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I. INTRODUCTION 



Black holes constitute one of the most exciting predictions of four-dimensional general rel- 
ativity (GR) and the astrophysical detections of them constitute a great advance in physics. 
Black holes further played a major role in establishing several mathematical aspects of GR 
and in understanding the physical interpretations of its solutions [H, 0, II]- For numerous 



gravity in higher dimensions has increased 
7J). The question of the existence of black 



reasons, the interest for general relativity and of 
considerably in the last years (see e.g. 01,0, 0, 
holes in higher dimensions therefore occurred naturally. One striking feature of the Einstein 
equations in more than four dimensions is that many uniqueness properties holding in four 
dimensions are lost. For example, black holes presenting an horizon of topology Sd-2 can 
be constructed explicitely in arbitrary dimensions [1]; their rotating generalisations U the 
Myers- Perry (MP) solutions, are characterized by several angular momenta. However several 
other types of solutions exist in higher dimensional gravity. Black strings are one of them; 
they denote string-like generalisations of 4-dimensional black holes in d-dimensional Einstein 
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gravity; in particular, they are characterized by an event horizon of topology Sa-3 x Si 
The simplest black strings are independant on the extra coordinate, say y. Later on it was 
realized that nonuniform solutions (depending on y) exist as well [ll|; these can further be 
done rotating (l2| . 



Among possible solutions in higher dimensional gravity, we should also mention black ring 



[13|, with toroidal horizon topology and black saturn 14(, where a black hole is surrounded 
by a black ring etc. These black object result from a balance between rotation along the 
torus and gravitational attraction and are known to be connected to the Myers Perry black 
hole in a limit where the rotation becomes high 13|, ll4| . However, it has been found recently 
that black saturns can also exist without rotation in Einstein Maxwell theory [IBj]. In this 
solution, it is the charge that balances gravity instead of rotation. 

As a consequence, the classification of localized objects (black holes, black strings or 
other) in higher dimensional gravity poses a serious mathematical problem. 

Recently, there has been a lot of interest for d-dimensional GR with a cosmological 
constant A. In the case of a negative cosmological constant, this interest is motivated by 
the correspondance between the gravitating fields in an AdS space-time and conformal field 



The propotype of higher-dimensional 
They generalize 

3 S3 



theory on the boundary of the AdS space-time [la . [17 
black holes with cosmological constant are known in an explicit form 
the solutions of [1,0]. Further exp licit solutions are known in non minimal model e.g. 
or in the presence of a dilaton 21, 22| . As far as we know, rotating black holes for the Maxwell 
field minimally coupled to Einstein gravity do not exist in a closed form and one has to rely 
on perturbative or numerical methods to construct them respectively in asymptotically flat 
' Q, anti-DeSitter [H, [2fJ and DeSitter [3 



space-times. The numerical construction 
251 ] uses the isotropic coordinates, making the comparaison with the explicit 
(usually expressed in Schwarzchild coordinate) not easy. In particular, the 




reported in 
solution of 

deformation of the MP analytic solution due to the electromagnetic field can hardly be 
estimated as explained later in the text. On the other hand, for A > 0, the occurence of a 
cosmological horizon leads to technical difficulties which do not allow a direct comparaison 
between the domain of existence of the MP solutions and those of charged black holes for 
both signs of the cosmological constant. 

The purpose of this paper is to reexamine the equations of rotating black holes in d- 
dimensional Einstein gravity with a cosmological constant and with a minimal coupling 
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to the electromagnetic field. We can then study the effects both of the sign of A and 
of the electromagnetic field on the MP solutions. The charged rotating solutions exist for 
angular momenta higher than some critical momenta. When this critical value is approached, 
the solution becomes extremal. However, we were able to address some extremal solution 
numerically by implementing a suitable set of boundary conditions. Extremal solution are 
interesting on their own, in particular, they are known to present symmetry enhancement 
[28| and play an important role in supergravity theory j29| . 

We consider odd values d, the ansatz for the metric is inspired from 0, [23[ . It contains 
(d + l)/2 killing vectors, namely time-translation and (d — l)/2 related to rotations in 
suitable two-planes of space. Assuming all angular momenta in the two-planes to be equal 
and completing the ansatz by an appropriate form for the electromagnetic fields, we can 
transform the Einstein-Maxwell into a system of six differential equations. The remaining 
arbitrariness in the coordinate system is fixed by using a Schwarzchild-like radial variable. 
Charged, rotating black holes can then be constructed numerically by fixing three parameters 
by hand, namely the values of the horizon, the angular velocity and the magnetic field at the 
horizon. The challenging question is to determine the domain of existence of the solutions 
in this three-parameter space. 

Although the equations are basically the same for both signs of the cosmological constant, 
the pattern of the solutions can be considerably different. 

The paper is organized as follow : in Sect. 2 we present the model, the ansatz and the 
relevant one- dimensional reduced effective action (the equations are written in an Appendix). 
In Sect 3 we discuss the explicit solutions known in some limits, namely the non rotating 
case and the vacuum case. The boundary conditions are presented in Sect. 4, together with 
some physical quantities characterizing the solutions. The numerical results are described 
in Sect. 5 and illustrated by several figures. In particular, the domain of existence of the 
black holes is emphasized for both A > and A < and compared to the vacuum case. 
Finally, some conclusions are drawn in Sect. 6. 



II. GENERAL FORMALISM 
A. The Action 

We consider the Einstein-Maxwell equations with a cosmological constant A in a 
d— dimensional spacetime. The equations follow from the variations of the action 

id 



d d xV^(R - 2A - F^Fn ~ TTT / d^xV-hK, (1) 

with respect to the metric and the electromagnetic field. Here Gd denotes the d-dimensional 
Newton constant and the units are chosen in such a way that Gd appears as an overal factor. 



The last term is the Gibbons-Hawking surface term 30]. It is required for the variational 
principle to be well-defined . The factor K represents the trace of the extrinsic curvature 
for the boundary DM. and h is the induced metric on the boundary. It is standard to 
parametrize the cosmological constant by means of an (anti)-DeSitter radius £ according to 
A = e(d — 2)(d — l)/(2£ 2 ); e = +1 and e = — 1 corresponds respectively to an asymptotic 
DeSitter and Anti-DeSitter space-time. 
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B. The ansatz 



To obtain rotating black hole solutions, representing charged U(l) generalizations of the 
MP solutions we consider space-times with odd dimensions: d = 2N + 1 and we parametrize 
the metric by means of an Ansatz previously used for asymptotically flat solutions [H, [23| 



ds 2 



-b{r)dt 2 



dr 2 



N-l /i-l 
i=l \j=0 



f(r, 

N /fc-1 



N /fc-1 



k=l \l=0 

N /fc-1 



fc=l \l=o 



yM 1 T c ° s2 ®i sin2 Okdfk 



_k=l \l=0 



(2) 



In the above formula 9q = and On = 7r/2 are assumed; the non trivial angles have 
9i G [0, 7r/2] for i = 1, . . . , N — 1, while ipk € [0, 2n] for k — 1, . . . , N. The corresponding 
space-time possesses N commuting Killing-vectors rjk = d Vk supplementing the standard 
invariance of time translations d t . The functions b, f, h, g,w depend on the variable r, the 
consistency of the ansatz imposes p(r) = g(r) — h(r). 

The most general Maxwell potential consistent with the symmetries of the line element 
(ED reads 



N /fc-1 



= V(r)dt + M 7 ") S ( ] I cos2 01 sin2 6kdlf>k 



(3) 



fc=l \l=o 



Where the electric and magnetic potentials V(r) and a v (r) depend on r. 

Substituting the ansatz , (J3J) into the action leads to an effective one-dimensional 
reduced Lagrangian 



S red J dr{Cf + 



(4) 



where the Einstein part C e J and the Maxwell part £ e J/ read respectively : 



Cf = (d- S)g^Jj((d - l)g -h) + \jjhbg^( b - + {d- 3) 9 -)(^ + (d - 3) 9 -) (5) 



Ud - 2)(d - D^Jhbg^g" + Ig^hMvfi 



(d-2)(d-l) ii-3) bh 
e- — -a 2 a — 

p 9 V / 



£$' = ^ (26 (2(d - 3)alh + fg 2 a' 2 ) - 2fg 2 h(wa' v + V') 2 ) ■ (6) 

The Einstein-Maxwell equations obtained from the ansatz ([S]),© can be obtained in a stan- 
dard way. For the numerical construction of the solutions, the radial variable r has to be 
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specified by fixing the "metric gauge"; we find it convenient to fix this arbitrariness by choos- 
ing g(r) = r 2 , defining a radial coordinate r of the Schwarzchild-type. The corresponding 
equations are written explicitely in an Appendix. Similarly to the asymptotically flat case 
23d . the electric potential can be eliminated from the equations by making use of the first 



integral (IBTj) . The cosmological constant can be arbitrarily redefined by means of a rescaling 
of the radial variable r and of the fields h, V and w. For the numerical analysis, we use this 
arbitrariness to choose r/, = 1 without loosing generality. 



III. KNOWN SOLUTIONS 

Solutions of the system of equations (}3"Tl) - (l3"oT) are known in two particular cases. 



(i) The de Sitter-Reissner-Nordstrom black holes [3l|, |32| are recovered in the limit w(r 

a<p(r) = : 

P ? \ 7 / \ r 2 2M q 2 

f(r) = 6(r) = l- c --_ + 



I 2 r d " 3 2{d - 2){d - 3)r 2 ( rf - 3 ) ' 
h(r) = g(r) = r 2 , V(r) = ^ _ ^ (7) 

where M and q are the mass and electric c harg e of the solution. 



(ii) The vacuum black holes discussed in [18J are recovered for a vanishing gauge fields 
(V = a v = 0) and 



r 



2 2ME 2Ma 2 , , s 9/J 2Ma 2 



where M and a are two constants related to the solution's mass and angular momentum and 
H = 1 + a 2 /£ 2 . These are generalisations of the Myers-Perry (MP) solutions [9[ to the case 
of non- vanishing cosmological constant. Here and in the following, Q = w(rh) will denote 
the angular velocity at the event horizon. 

Some properties of the MP solutions are useful for the discussion of charged rotating black 
holes in the case A > 0. If we consider uncharged black holes solutions with fixed and r c , 
the explicit form provided by the Myers-Perry solutions reveals that rotating solutions exist 
up to a maximal value of Q. For d = 5, we find 



/ 2 _ r cr\{r 2 c + r 2 ) I 1 

1 ' r 2 + 2r 2 r 2 rl + r c 2 + 2r 2 ' & r 2 + 2r\ K ) 

For Q — > Vt max the event horizon becomes a double root of the function /, i.e. the horizon 
becomes extremal. In fact for Q > Vt max , a third horizon occurs in the interval [r^, r c ] in the 
neighbourhood of r^. The values fl ma x and the corresponding value of the DeSitter radius 
l/£ 2 are reported as functions of the cosmological horizon r c for three different values of 
on FigfTJ 

For fixed values of r^, r c , the value of £ 2 varies only a little for Q, e [0, n maa; ] ; e.g. for 
Th — l,r c — 3 we have 1/10 < l/£ 2 < 1/11. The mass E and the angular momentum J of 



5 



1.0 



0.6 



0.4 



0.2 



0.0 



I 
I 

\ 

\ 




\ 

\ 






dash : r h = 0.5 


A 

\ 

v \ 

\ \ 
v n 


solid : r h = 1 .0 
dots : r h = 2.0 







■ i i i 



10 



FIG. 1: The value O maa; and the parameter l/£ 2 as functions of r c for the MP solution and for 
r h = 0.5,1.0,2.0 



the Myers- Perry solutions are given by [33j, [34 



E 



Vd-o „ r ,d 



-2 a 
2 e 2£ 2 



v: 



d-2 



4vrG 2 



-Ma 



(10) 



where Vd-i denotes the area of the d — 2 dimensional sphere. Throughout this paper, we 
will use the expression E above to characterize the mass of the solutions although a more 
elaborated formalism exists 35|, [36] . Once applied to the MP solutions 38| , it leads to an 
expression of the form = E^ — eE where E^ is interpreted as the Casimir energy of 
the dual Conformal Field Theory living on the d— 1 -dimensional boundary of the space-time 
The contribution Ec appears only for odd values of d and depends on the space- 
time dimension and of the cosmological constant 37]; for d — 5, one finds E^ = -^f^ 



IV. BOUNDARY CONDITIONS AND PHYSICAL QUANTITIES 
A. Constraints of regularity about the horizon 

We are interested in black hole solutions, with an horizon located at r = rh, i.e. with 
b{rh) = fi^h) = 0. Expanding the different functions in the neighbourhood of the horizon 

b(r) = 6i(r - r h ) + 0(r - r h ) 2 , f(r) = f x (r - r h ) + 0(r - r h ) 2 , 
h(r) = ho + hx(r - r h ) + 0(r - r h ) 2 , w(r) = w h + w x (r - r h ) + 0(r - r h ) 2 , 
a ip (r) = a + a 1 (r-r h ) + O(r-r h ) 2 , V(r) = V + V x {r - r h ) + 0(r - r h ) 2 (11) 



(3 



and inserting the expansions (fTTj) into the equations (|3T]) - (!36|) leads to terms proportional to 
l/{r — rh}. For the solutions to be regular at the event horizon the coefficients of all these 
singular terms should vanish. Since our numerical analysis will deal mainly with the case 
d = 5, e = 1, we write these rather cumbersome conditions in this case only. The equation 
for the function h leads to the condition Ti(r = r^) = with 

r x (r) = 8b'h 2 (12al + 7h) + 4rb'hh'(12al + 5h) 



- 32&'/iV + 8r s b'h{f'h - Ah') 
24 

+ 1hr\— b'h 2 - f'{A(a') 2 hw 2 + 8a'hwV + b'h! + 5h 2 (w') 2 + Ah(V') 2 )) 



24 

+ h'r 5 (-^b'h 2 - f'(4(a' v ) 2 hw 2 + 8a' v hwV - b'h' - h 2 (w') 2 + Ah(V') 2 )) (12) 

where, for brevity, we dropped the dependence on r in the various functions. The value 
f'{rh) appearing in the condition ri(r^) = can be extracted from the equation for /, 
giving 

f(r x _ 4hh 6rj/i 2 + 8r 2 - 12a 2 - 5h 

J[h) ~ rl ^b 1 h + r h {Ah Q {V l + a l w,f-b 1 h l -h 2 wl) 1 ) 

In the same way, the two Maxwell equations lead to a single condition T 2 (r = 77J = with 

T 2 (r) ee Aa v b'h + r A f\hw'V + a'^hww' - a'^b') (14) 

Similar conditions should be imposed at the cosmological horizon r c as well : r x (r c ) = 

r 2 (r c ) = 0. 



B. Physical quantities at the horizon 

The Hawking temperature and the event horizon area of these solutions can be obtained 
in a standard way, leading to 

Th = ^, A H = V d . 2 rt 2 . (15) 

The mass and the angular velocity at the horizon are also usefull. They are defined by 
means of the appropriate Komar integrals over the sphere Sjj at the event horizon. It leads 
to 

Vd-2 I fhg 2 , , Vd-2 n I fg 2 h 3 ,. 

These quantities can be evaluated from the numerical solutions. Similar quantities, say 
Mc, Jc can be associated with the cosmological horizon. In fact, the masses and angular 
velocities at the two horizons can be related to each other by formulas of the Smarr-type. 
Considering a Killing vector K^, integrating the two sides of the identity 

VV^ = -WjCn = -8vr(Tt - J-JTSfiKp (17) 

over the (truncated) hyper- volume E covering the space-like region between the two horizons 
and using both Einstein equations and Stokes theorem appropriately result into the identities 
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we are interested in. Choosing first K = one gets 



Jh-V Cj 



d-2 



1W 



\r=r h , 



(18) 



expressing the conservation of the total angular momentum. The identity associated with 
the time-translation invariance (K = dt) is a bit more involved, it leads to 



M c + §cQc -M H - $ H Q H 



d-2 



fhg< 



dr 



(19) 



where Qh (the electric charge at the horizon 39J) and (the electrostatic potential) are 
defined respectively by 



Q 



\r=r h j 



H 



d 



-Qa 



ipj\r=r h 



(20) 



and similarly for Qr, Similarly to 

ity k 

d-3 



the mass at the horizon can further be expressed 
in terms of the surface gravity k and angular momentum according to 



-M 



u 



1 a d ~ 1 ^T 

-kAh H \IJh 



d-2 8xG d "~~" ' 2 

The various identities above were checked on our numerical solutions. 



(21) 



C. Asymptotics and global charges 

The asymptotic form of the solutions to the Einstein equations can be constructed after 
some algebraic manipulations, leading to 

b(r) = -e£ + 1 + -^3 + 0(l/r™-% f(r) = + 1 + ^ + 0(l/r^), (22) 



h(r) = r 2 (l + + 0(l/r 2d ~ 4 )), w(r) = 4=1 + 0(l/r M " 4 ). 



They depend on three arbitrary constants a, (3 and J. The choice of the gauge metric g = r 2 
in the ansatz ([2]) and the resulting functions /, b, h have the advantage to present a direct 
connection with the closed form of the MP rotating solution. The asymptotic expression of 
the gauge potential is similar to the asymptotically flat case 

V(r) = + 0(l/r 2 ^), a v {r) = + 0(l/r 2 ^). (23) 

Although it is not clear if the thermodynamical properties are well defined in the presence 
of the positive cosmological horizon, one may still define them in the standard way. In this 
section, we follow the lines of j25| to present the global charges characterizing the solutions 
asymptotically. [35|, 36[. The mass-energy E of the solutions and the angular momentum J 



associated with an angular direction are respectively given by 
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The above relations can be obtained namely by using a background subtraction approach or 
the counterterm formalism 35|, 136J1. Here do not take into account the contribution of the 
Casimir energy of the solution 

The electric charge and the magnetic moment of the solutions are given by 




Q 



(d-3)V d -2 
47rG n 



-q, n 



(d-3)V c 
AnG n 



d-2 



(25) 



D. Boundary value problem 

1. Case A < 

In the case of asymtotically ADS black holes, the formulation of the boundary problem is 
simple. Fixing A < 0, the equations can be integrated in one step on the interval r G [r^, oo] 
by imposing the boundary conditions 

/ = ,6 = 0,«; = «; h ,y = 0,a; = a;,r 1 = 0,r 2 = for r = r h , (26) 

,2 

b= r — + l + 0(l/r 2 ) , h = r 2 , (r 2 a<j> y = 0, (r A w)' = for r -> oo . (27) 

The arbitrary additive constant of the electric potential has been used to set V(rh) = 0. 
The parameters Wh, a' h have to be fixed by hand (alternatively, a'^ = a' h can be replaced by 
a <p = o/i)- The asymptotic decay of the electric potential V and of the metric function / can 
then be used for crosschecking the numerical solutions. 

Extremal black holes are characterized by a double zero of the functions f,b at r = rv In 
the case of de Sitter- Reissner-Nordstrom black holes, analytical form follows from (171), the 
parameters M, q are then fixed as functions of r/j. In the case of rotating solutions, the 
regularity conditions at r = r^, leads to a set of rather involved conditions at r = rh- 



2. Case A > 

The positive cosmological constant leads to the occurence of a A-depending cosmological 
horizon at r = r c with < r c < oo and where f(r c ) = b(r c ) = 0. This creates several 
difficulties, namely : (i) the point r c constitutes a singular point of the equations, (i) the point 
r c is not known a priori as a function of A, (iii) imposing regularity at both horizons and 
r c needs twelfe boundary conditions, this number exceeds the eleven boundary conditions 
allowed by orders of the six equations (for instance, one of first order and five of second 
order) . 

In order to overcome these difficulties, we supplement the system (T5T|) . (1561) by the trivial 
equation d£ 2 /dx = 0. The new system then needs twelve boundary conditions to be well 
posed. The strategy consists in solving the equations in two steps. First we choose by hand 
the event horizon and the cosmological horizon r c and integrate the equations on the 
interval [r^,r c ] by imposing the regularity conditions at the two horizons. It turns out that 
this needs exactly twelve boundary conditions, namely : 

/ = , 6 = 0, b' = 1 , w = w h , V = , a'=a , h ,V 1 = 0, T 2 = for r = r h (28) 
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/ = , b = , Ti = , T 2 = for r = r c (29) 

The functions r 12 are defined above. Here, we choose the arbitrary rescaling of the time 
variable t in order to impose b'{rh) = 1- The function b(r) is therefore not normalized 
according to ff22l) . the appropriate normalisation, say b — > b/M 2 , can be obtained only 
after the second step. The constants Wh and a! h are arbitrary and control the total angular 
momentum and the (electric and magnetic) charges of the black hole. The numerical value 
of I 2 is determined by the first step, together with the values of all the different fields at 
r = r c . 

The second step consists in integrating the equations on [r c , oo]. The value of I 2 is known 
from the first step and the numerical values of the fields at r = r c can be used as a set of 
Cauchy data to solve the equations from r = r c to r > r c . 

This method presents several features of inconvenience. Namely: (i) a systematic analysis 
of the solution for a fixed value of the cosmological constant cannot be performed. Fortu- 
nately, the numerical value of £ 2 depends only a little on w{r h ) and a v (r h ) once r h ,r c are 
fixed, (ii) The appropriate normalisation of b,w,V leading to an asymptotically DeSitter 
metric is only known a posteriori. 

h b w v V 

A/; 2 A4 M b 

In particular it turns out impossible to study solutions with a fixed value of the angular 
velocity at the horizon Q (since Q = w(rh)/Nb) and varying the horizon. In the numerical 
analysis, we put the emphasis on the families of solutions obtained when varying Q for fixed 
values of the horizons r^, r c and of the parameter a' h . 



V. NUMERICAL RESULTS 

The system of equations (|HT|) - (|5S1) does not admit, to our knowledge, explicit solutions 
for generic values of Wh and a' h . We therefore rely on a numerical method to construct 
solutions. We solved the equations in the case d — 5 and we hope that this case catches 
the qualitative properties of the pattern of the solutions available for d > 5. The numerical 
solver Colsys [40j was used to obtain the results. 



A. Charged solutions for A < 

Charged 5-dimensional rotating black hole solutions have been presented in details in 
251 ] where the coordinate freedom is fixed by chosing the isotropic coordinate, say y. The 
solutions are constructed as function of y and their pattern is discussed mainly for fixed 
horizon at y = y^. With this definition of the radial coordinate, the pattern of the solutions 
available for fixed yh present several branches. For instance, solutions with different masses 
but same yh and Q can be found. Unfortunately, the solutions can hardly be compared with 
the the explicit form (jSJ) of the Myers-Perry solutions which hold with the Schwarzschild 
(or the Boyer-Linquist) coordinate. For this reason, and for the purpose of comparison of 
the solutions between the Anti-DeSitter and DeSitter cases, we present some solutions with 
A < in the Schwarzschild coordinate g = r 2 and determine the pattern of solutions available 
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FIG. 2: The Mass E, the angular momentum J and the value f\ as functions of 0, for the Myers- 
Perry solution (solid) and the charged solution with A' = 0.5 (dashed) for = 1 and I 2 = 0.1. 
The Hawking temperature Th is represented by the dotted lines 

with Th fixed. For all this section we set r^ — l and I 2 = 0.1. Several physical parameters 
characterizing the charged rotating black holes with r/j = 1, a' h = 0.5 are represented on 
Fig. [2] by the dashed lines. The corresponding parameters associated with the MP solutions 
are represented on this figure by the solid lines. As already pointed out, vacuum (or MP) 
solutions exist for positive values of Q up to a maximal value fi max , with the choice of 
the parameters of the figure, we have VL max ~ 0.78. Our numerical integration of the full 
equations reveals that charged black holes exist only for a finite interval of the horizon 
angular velocity : Q G [Q a , &>b\, with Q a > and < Vl max . The values Vt a ^ depend on the 
strength of the electromagnetic fields. In the case of Fig. [2] we find Q a ~ 0.22, f2& ~ 0.765. 
In the limits Q — > Q a ,b, extremal solutions are approached, i.e. corresponding to f\ = b\ = 0. 
As a consequence, the Hawking temperature Th also tends to zero in these extremal limits. 
It is represented by the dotted lines on Fig. [2j Several parameters characterizing the 
electromagnetic fields of the solutions have also been computed numerically. The electric 
charge Q and the magnetic moment fi, vary slowly as functions of the horizon velocity Q, as 
shown by Fig. [31 Note that that the numerical values of Q and V(oo) are very close to each 
other and similarly for fi and Ah- Profiles of a typical ADS rotating black hole is presented 
on Fig. |H We managed to construct rotating extremal black holes numerically for small 
values of Q. The profile of such a solution is presented in Fig. [5] corresponding to Q = 0.09. 
The numerical integration is delicate; a systematic construction if extremal solutions is not 
aimed in this paper but would require an appropriate reparametrisation of the metric. 
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FIG. 3: Parameters of the electromagnetic field of charged ADS-BH as function of £1 for rh = 1, 
ai = 0.5. 



B. Remarks on the Isotropic coordinate 

We would like to stress the main differences between the patterns of the solutions once 
constructed by using the isotropic coordinate and the Schwarzschild coordinate. We make 
these remarks in the vacuum, asymptotically flat case : A = 0. In the present paper, we use 
the Schwarschild coordinate r, denoting the angular velocity at the horizon as Q = w(rh) 
fixing the horizon r^, we have solutions with a regular horizon for fl < ^— All physical 
quantities, the mass in particular are monotonic functions of Q. In contrast, the isotropic 
coordinate y is used e.g. in j25[, we have 

y = expK(r) , K'{r) = (30) 

rvf 

While families of solutions (MP) take a rather simple form in the Schwarzschild coordinate 
it is not the case with the isotropic coordinate. Considering the family of solutions with 
fixed horizon y^ in the isotropic coodinate, leads to a rather involved pattern in the r^, Q 
plane. Namely it can be that two solutions with equal and fl can exist but with different 
masses. This is illustrated on Fig. [6] where two solutions corresponding to « 1 and pa 2 
and fl = 0.33 have both y^ = 0.5 once expressed in the isotropic coordinate. 



12 




FIG. 6: The profiles of two black hole solutions corresponding to = 0.6, 0, = 0.33 



C. Charged solutions for A > 

For A > 0, the presence of a Maxwell field changes the pattern of solutions corresponding 
to the vacuum (or MP) solutions. Along with the case A < 0, our numerical analysis 
strongly suggests that charged-rotating black holes exist only for a finite interval of the 
horizon angular velocity Q, i.e. for Q £ [fl OJ f]j,] with < Q a < Qb < Q max . Here Vt max 
denotes the maximal horizon velocity for the MP solution and the values Q a ^ depend of the 
strength of the electromagnetic field. It must be stressed that the numerical integration 
of the equations become extremely tedious in the region of the parameters which allow to 
approach the lower limit Q ~ Q a . The reason of this difficulty comes from the fact that, when 
the shooting parameter Wh decreases, the branch of solutions stops as some minimal value, 
say Wh,min- Then, a second branch of solution emerge for Wh > Wh,min- The two branches 
coincide it the limit Wh Wh.min- Following the second branch shows that it terminate into 
a critical value of Wh, c and fx — > for Wh — > Wh,c- This is illustrated on Fig. [7] where a few 
parameters that are substantially different on the two branches are represented as functions 
of Wh- In all figures of this section, the values r h = l,r c = 3 are used corresponding to 
£ 2 ~ 10. In fact the second branch is represented for a' h = 0.5 (where we find Wh t min ~ 0.459 
and Wh, c ~ 0.67). A second branch exist also for the other lines presented on the figure (i.e. 
corresponding to a' h = 0.1) but it turned out even more difficult to construct it with a good 
accuracy than for the case a' h = 0.5 and we refrain to representing it. On the other side on 
the graphic, solutions seem to exist for arbitrarily large values of the parameter Wh- 

So, once described in terms of Wh, the pattern of solution looks complicated. However, 
it is more relevant to plot the different quantities as functions of the physically meaningfull 
parameter Q representing the angular velocity at the horizon. This is obtained by renormal- 
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12 3 4 



w h 

FIG. 7: The values f\, V\ and Vinf = V(oo) are reported as functions of Wh for a' h = 0,0.1 and 
a' h = 0.5 

izing the fields b, w, V appropriately so that the metric become asymptocally DeSitter. 

Several parameters characterizing the metric functions at the horizon are reported on Fig. 
[HI as function of Q. The results summarized on this graphic suggest that the solutions exist 
for Q e [Q a , fife], and that in the two limits Q —>■ an extremal black hole is approached 
at the event horizon r^. Physical parameters like the energy E and the angular momentum 
J, which are related to the asymptotic charges, are reported on Fig. [91 The two figures 
[U M further show that, for a' h fixed, the parameters of the charged black hole deviate more 
substantially from the corresponding Myers-Perry solutions in the region of small values 
of Q. We note that, for non-vanishing electromagnetic field there is only one solution for 
given Q, r^. This contrasts with the description in term the isotropic coodinate where in 
general two solutions can be associated to one value of Q (see Fig.l of (25|). Finally several 
parameters related to the Maxwell fields V and a v are reported of FigJTDl We see that 
the values of fi,A which hardly differ for small a' h are lifted when the amplitude of the 
electromagnetic field increases (and similarly for Q, Kn/)- 

Profiles of two typical solutions having a' h = 0.5 are represented on Figs. [TTyT2l repectively 
for Q = 0.62, i.e. close to the limit flj, ~ 0.64 and for an intermediate value Q = 0.44 
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FIG. 8: The value f\ and go are reported as functions of for a' h = 0, 0.1 and a' h = 0.5 



VI. CONCLUSIONS 

In this paper, numerical aguments are given for the existence of five-dimensional charged 
and rotating black holes with a positive cosmological constant. Two main features of these 
solutions are that (i) they are characterized by at least two horizons, an event horizon and 
a cosmological horizon r c ; for fixed r^,r c , we find a families of black holes characterized by the 
parameters Q and a' h determining respectively the angular velocity and the magnetic fields 
at the event horizon. We presented several arguments indicating that black hole solutions 
exist only for Q taking values in a finite interval, say [fl a , fib]. The limiting values depend on 
the strength of the electromagnetic field (in our case this is controlled by parameter a' h ) and 
determine the domain of existence of the black holes in the fl, a! h plane (with r^, r c fixed). 
The solutions end up into extremal black holes when the values fi 0j & are approached. In 
the limit a' h —* 0, the MP black holes are recovered, existing for Q £ [0, f2 maa; ] with fl max 
known analytically. This feature is common to both signs of the cosmological constant once 
solution are constructed with a Schwarschild -like coordinate. Smarr-type formulas have 
been obtained which relate some physical quantitities associated with the event horizon and 
the cosmological horizon in the case of a positive cosmological constant. 
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FIG. 9: The Energy of the solutions and the angular momentum J are reported as functions of Q 
for a' h = 0;0.1 and a' h = 0.5 

VII. APPENDIX: THE EQUATIONS 

In this Appendix, we present the equations associated to the Einstein-Maxwell La- 
grangian. First, the equations determining the metric : 

fl , f ( rh ,2 . 2r T/ / 2 . 4 ™ / T w . h' rb' 1 w 2 /2 6' 

/ + (Z^y^-^ +T r + ^ a ^ + X (1 -26 } - 2r( 7I-T )a - + 6 (31) 

( "- 1)(rf - 4) ^ ' 1 :((3rf-5)/, + 4( C / + l)a^-( C /-l)V)+ (d - 1)r -- 



6" + — ?— f 4(5 - 2d)wa'V + ^-^-b'h' - 2( - d ~ h ' - 2{2d - 5)V' 2 (32) 
a — 2 \ ^ 2h 2rh 

+1(3 - 2d) W 2 - 2(£ + (2d - + (d - 2) - Q + 

(d-3)b, „ (d-3)6/d-l , \ (d-l)(d-2)6\ 

; -(12a 2 +/i) + ^ 5^- __ + 4-d +^ ^ '-) =0, 
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FIG. 10: Electromagnetic parameters as functions of 0, for a' h 0.1 and a' h = 0.5 

1 ( {2d-h)h\ j2 2h Ahw _ (d-2)h\ h< _ f 
+ (d-2)\ 2b + b + b a * V 2 [ h f } 

(d-3) (rf-3) 2 ^ o ,W 2 o ^ 2 (d-3)fr (rf- 3 )(2rf-3)/z 2 

* + ^ + 2{ — + 2d ~ 5)a - ~ ~^b— b ^7 

12(rf-3)a 2 /i (rf_3)(rf_4)/i (d-l)hd-2 d-3 



2 



a ^ ft r 2 \ b j n J 

and, for the Maxwell fields 

T/" W h' > a. ™ ' -L 1 /' 

^ -^ a , + I a ^ + 2 ( ^ -b + J + J )V 

+ (1 + ^)a>' + + 2(rf -^ W = 0, 



„ l,2(d-3) b' f h\ , h. , , 2(d-3)a ip h 

It can easily be seen that the equations of motion have the first integral 



g^^(wa' v + V') = (d-3) q . 
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A' = 0.5 , n = 0.62 




1 2 3 456789 10 

r 

FIG. 11: The profile of the metric and Maxwell function for r c = 3,r/j = 1 for a' h = 0.5 and 
O = 0.62 




1 2 3 456789 10 
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FIG. 12: The profile of the metric and Maxwell function for r c = 3, = 1 for a' h = 0.5 and 
n = 0.44 
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